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Nowa metoda rozkładu liczb na czynniki 
pierwsze. 


Wszelkie znane dotąd w matematyce sposoby rozkładu wielkich liczb 
na czynniki pierwsze nie są właściwie metodami rozkładu, lecz metodami 
wynajdywania form kwadratowych i linjowych dla dzielników 
poszukiwanych. Bardzo dowcipne i subtelne sposoby matematyków tej miary, 
co Lagrange, Legendre, Lucas, a zwłaszcza Czebyszew ustalają przedewszyst- 
kim ciąg form kwadratowych, jedynie możliwych dla czynników liczby bada- 
nej, — a zarazem ciąg form linjowych, odpowiadających każdej kwadratowej. 
Kombinując następnie po dwie formy kwadratowe w różnych zestawieniach, 
wyłączamy z pomiędzy odpowiednich im linjowych te, które nie należą jed- 
nocześnie do obu kwadratowych, zestawionych przez nas. W ten sposób otrzy- 
mujemy formę linjową coraz ogólniejszą, w której * spółczynniki liczbowe 
wciąż wzrastają. Zatrzymujemy się wreszcie na takiej, która da nam możli- 
wie najmniejszą ilość liczb pierwszych jej odpowiadających, i w końcu wy- 
konywamy po kolei dzielenie próbne liczby badanej przez wszystkie 
znalezione przez nas możliwe jej czynniki. 

Postępowanie takie dalekie jest od doskonałości nie tylko dlatego, że 
nie daje nam czynników gotowych, że je wyczuwa zaledwie. Większą 
jeszcze jego wadą jest to, że zastosowywane tu działania skompliko- 
wane i długie, a wymagające uwagi skupionej, nie podlegają żadnemu : 
sprawdzaniu się w biegu samej roboty. Możliwe więc są omyłki przy- 
padkowe, unicestwiające wynik ostateczny. Dość zwrócić uwagę na wzór za- 
sadniczy Czebyszewa, na podstawie którego otrzymujemy szereg form kwa- 
dratowych na dzielniki liczby badanej. Jest on tak skomplikowany,*że wszel- 
kie działania stąd wypływające należy przerobić powtórnie, lub dać do prze- 
robienia osobie innej, by mieć pewność, że punkt wyjścia nie zawiera 
omyłek, 

To, eo podaję poniżej, nazwałem metodą rozkładu liezb na czyn- 
niki, gdyż mając badaną liczbę N, wykonywam pewne działania algiebraiczne 
i arytmetyczne, które prowadzą bez żadnego dzielenia próbnego 
albo do otrzymania gotowych dzielników liczbowych, albo do pewności, 
że liczba N jest pierwszą. Prócz tego, jeśli N okaże się złożoną, mogę 
orzec z łatwością i a priori, wiele mianowicie zawiera ona różnych czyn- 
ników pierwszych. Omyłki przypadkowe wyłączone są absolutnie 
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w tej metodzie, gdyż mamy tu działania najprostsze i sprawdzenie natych- 
miastowe roboty głównej, wyznaczającej liczbę dzielników. Wreszcie, jeśli 
chodzi o szybkość, to np. liczbę 8520191, eytowaną i rozkładaną przez 
Czebyszewa, zbadałem przy pomocy mojej metody w ciągu 1!/, godziny, po- 
siłkując się zaś arytmometrem— wciągu 40 minut. Tymczasem, przerobienie 
tego wszystkiego co wykonywa sam Czebyszew *), zajęło mi przeszło 12 
godzin, i prócz tego — robota cała wymaga przerobienia powtórnego. 

Wskazuję tu na wyższość mojej metody; słusznym jest zwrócić uwagę 
i na jej słabą stronę, a mianowicie, że, jak dotąd, może być ona stosowana 
praktycznie tylko do liczb kształtu 10m--1. Co zaś do kształtu 10m--3, to 
jakkolwiek podaję teoretyczne rozwiązanie kwestji**), lecz wymaga ono 
przejścia do kwadratu liczby N, co prowadzi do liczb zbyt wielkich i pozba- 
wia metodę jednej z głównych jej zalet, t: j. szybkości. Dla liczb kształtu 
10m-3 należy więc szukać metody innej, krótszej. Zdecydowałem się jednak 
na publikację wyników już osiągniętych, w nadziei zachęcenia matematyków 
do opracowania i udoskonalenia dalszego metody tu podanej, jako skierowu- 
jącej kwestję omawianą na właściwe tory. 

Na wynalezienie jej naprowadziły mię badania nad włówności Gi 
uogólnionych ciągów Fibonaccie'go, ogłoszone już drukiem ***). 

Ciekawych niezmiernie badań ciągu Fibonacciego dokonał E. Lucas 
(Bolletino Boncompagni. Tom X. Str. 129. Sur diverses questions d'arithmć- 
tique supérieure), gdzie wspomina nawiasowo o uogólnionym ciągu Fibonac- 
ciego i jego określniku, nie robiąc jednak żadnego stąd wniosku o możliwo- 
ści zastosowania go do rozkładu liczb na czynniki pierwsze. Lucas rozwinął 
następnie swą analizę, stosując ją do trzech typów ciągów: Fibonacciego, 
Fermata i Pella (Théorie des fonctions numériques simplement periodiques. 
American Journal of Mathematics. 1878. Tom I. Str. 184). 

Przypomnę tu w streszezeniu wyniki główne badań moich, potrzebne 
do zrozumienia metody. przedstawionej poniżej. 

Cheąc, by ciąg uogólniony nie był skracalny, SE dwie wielkości 
dodatnie: aib względnie pierwsze, i otrzymujemy postać algiebraicz- 
ną takich szeregów 


a, b, a+b, a--20, 2a--3b, 3a+-50, 5a--86, ... (1) 


gdzie, oczywiście, spółczynniki przy ai b w każdym wyrazie poszczegól- 
nym są kolejnemi wyrazami ciągu Fibonacci'ego 


Ime NEPA AE AEO CENY ZCZEC EE (2) 


Ciąg liczbowy, czy też rozważany w kształcie algiebraicznym, należy 
rozumieć zawsze, jako doprowadzony do właściwego swego początku, co 
jest uwarunkowane przez nierówność 


*) Czebyszew. Sobranje Soczinienij. Tom I. Str. 73. 


2! **) R. Niewiadomski. Nowyj metod razłożenja czisieł na pierwonaczalnyje 
mnożiteli. Moskwa 1912. 


**) ` Wiadomości Matematyczne, Tom XV. 1911. "Sty, 230. 
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b—a>a czyli v>2a (3) 


Przy zachowaniu tego warunku ciąg rozpoczynać się będzie od liczby mo- 
żliwie najmniejszej. 

Dla każdego ciągu uogólnionego istnieje pewna liczba stała, niepa- 
rzysta D, czyli określnik szeregu, która równa się 


DEE UŻy— Uy1 . Uy41) (4) 


Przy N parzystym otrzymujemy przed nawiasem znak --, przy N zaś 
nieparzystym znak —-. 
Kształt algiebraiczny określnika ciągu typowego (1) będzie 


D=v*—a(a--b) =b*—ab—a* (5) 


czyli odpowiada mu pewna specjalna forma kwadratowa. 

Określnik D jest liczbą pierwszą względem wszystkich 
wyrazów szeregu, lub szeregów mu odpowiadających (o ile a jest pierw- 
sze względem b). 

Jeśli mamy jakikolwiek ciąg (1) o określniku D, to istnieje co naj- 
mniej drugi jeszcze ciąg o tymże określniku, a mianowicie 


b—2a, 2b— 3a, 3b—5a, 5b—8a, 8b—13a, 13b—2la, ... (6) 


Oba te ciągi (1) i (6) mają jednakowy przedwyraz: U,=b—a 
będziemy je nazywali sprzężonemi. Z dwuch sprzężonych ciągów nazy- 
wać będziemy zasadniczym ten, którego wyrazy początkowe są względnie 
mniejsze. 

Dwa jedynie ciągi nie posiadają sprzężonych: 1) ciąg Fibonaeci'ego 
(e gdzie Us=l, 1 DI=]; 2).6iąg 1, 3, 4%, 7, 11,518, ... w. którym 
U5=2 i D*=5. 

Określniki D mogą być tylko kształtu 10ml, o ile nie zawie- 
rają czynnika 5, który zresztą występować może w określniku jedynie w po- 
tędze pierwszej. 

Badając kształt ogólny ciągu sprzężonego (6). widzimy, iż spółezyn- 
niki przy a i b są kolejnemi wyrazami ciągu (2), poczynając od wyrazu trze- 
ciego. Rozważywszy zasadę formowania ciągu (6) z ciągu (1), znalazłem, iż 
zupełnie identycznie powstaje ciąg o określniku złożonym D=dyd, z dwuch 
innych ciągów, których określniki równają się d, i d,. Mianowicie, jeśli 
pierwsze wyrazy ciągu mnożonego są a i b, to postać ogólna wyrazów ciągu 
złożonego będzie 
sę pb—qa (7) 


gdzie spółczynniki zmienne p i q są wyrazami kolejnemi ciągu mnożącego. 
Naturalnie, role ciągów mnożonego i mnożącego mogą się nawzajem zamie- 
niać. W przypadku szezególnym, ciągi (6) i (1) mają ten sam określnik, po- 
nieważ określnik ciągu mnożącego (Fibonacciego) równa się 1. 
Przypuśćmy, że dane ciągi o określnikach d, i d, mają postać 
(o) a, b, a+b, a+2b, 24+3b, ... gdzie d,=b—ab—a* 


(p) P, 9, P+9, P+-29, 2p--3q, ... gdzie dy =q*—pq—p*.. 
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Kształt ciągu o określniku złożonym D=dyd, powinien być 
"(0 pb—qa,. qb—(p--9)a; (P--9)b—(0--29)a, (p--24)6—(2p+-39)a, ... 


Ciąg (y) jest, oczywiście, jednym z uogólnionych ciągów Fibonacei'ego, 
ponieważ spółczynniki przy b i a w każdym jego wyrazie również tworzą 
taki ciąg, w którym suma dwuch wyrazów kolejnych równa się trzeciemu. 

Określnik ciągu (y) jest 


D„=|q0—(--9)a]*—(pl—-q0))(p--06—(0-P2g)aj (8) 
. 5 
Wystarcza sprawdzić, że zachodzi tu tożsamość z iloczynem 


(b6—ab—a*)(q*—pq —p*) (8a) 


Zachodzi to w rzeczy samej, jak łatwo się przekonać po wykonaniu 
mnożeń odpowiednich. 

A więc, iloczyn form kwadratowych kształtu (5) b—ab—a* 
ma taką samą formę kwadratową (5), zupełnie analogicznie 
do znanej własności form a*--b. 

Na rozwiązaniu odwrotnego zagadnienia drogą algie- 
braiezną opiera się moja metoda rozkładu liczb określni- 
kowych na czynniki pierwsze. 

Ponieważ przekonaliśmy się, że istnieje równość (8)=(8a), możemy 
więc wnioskować w związku z tym, co dowiedzione zostało powyżej, iż określ- 
nikami ciągów rozważanych mogą być jedynie: 1) liczba 5; 2) wszystkie 
liczby pierwsze kształtu 10m+1; 3) wszystkie kombinacje czynników grupy 
poprzedniej, z powtórzeniami lub bez takowych, i w kombinacji z liczbą 
5 w potędze pierwszej. 

Tablicę określników, doprowadzoną prawie do 1000, podałem w innym 
miejseu *). Tamże znajdzie czytelnik dowód tego, co następuje: 

Szereg mnożący o określniku 5 nie powięksża liczby ciągów złożonego, 
jak również nie ma wpływu na ten rezultat potęga czynników, zawarty ch 
w danym złożonym. Wogóle, dwom różnym określnikom 


o a A 
D=5. s% P, ul, wu, i D; =s. t, u. 0... 


odpowiada ta sama liezba ciągów, sprzężonych parami **). 

Dowodzę też sposobem popularnym ob że jeśli dany określnik złożony 
zawiera m różnych czynników pierwszych, nie rachując czynnika 5, to wszyst- 
kich ciągów o tym samym określniku mamy 2”, czyli 2%-! par ciągów sprzę- 
żonych. Ponieważ jednak wyrażenie (5) jest prócz tego formą kwadratową, 
zatym rezultat znaleziony przezemnie wypływa bezpośrednio z teorji tych 
form, gdzie dowiedziono, iż, mając N=a.f. 4... zawierające m czynników 
pierwszych, możemy N przedstawić w postaci 2% różnych form kwadrato- 


*) Wiadomości Matematyczne. Tom XV. 1911. Str. 239. 
**) Tamże. Str. 243—246. 


***) R. Niewiadomski. Nowyj Metod rozłożenja czisieł na pierwonaczalnyje 
mnożiteli. Moskwa 1912. . ; i .) 
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wych, czyli otrzymamy 2% różnych pierwiastków na a i 5, t.j: 2” różnych 
ciągów typu Fibonacciego. 

Gdy D jest liczbą złożoną, niezbędnym jest mieć przynajmniej dwie 
pary ciągów mu odpowiadających, by zastosować moją metodę rozkładu. 

Przypomnę tu, co pokazałem w innym miejscu na przykładzie liczbo- 
wym *), w jaki mianowicie sposób z dwuch par ciągów o określnikach d, i d, otrzy- 
muje się dwie pary nowych ciągów, których określnik równa się dd, Weź- 
my d,=11 i d,=19. 


sk zasadniczy: «1, 4; 5, 9; 14, 23, 37) (a) 
ża sprzężony: *«2, 2: 5500773. 12,: 19; 31, 550;04- «(B) 


AR zasadniczy: - 1, 5, 6; 11, 17, 28, 45, .. (1) 
de: sprzężony: 3, 7, 10, 17, 27, 44, 71, a () 
U —13 zasadniczy: 5,.18, 23, 41, 64, 105, 169, ... (s) 
SĘ | 

D=209; | sprzężony: 8, 21, 29; 50, 79, 129, -208, ... (6) 
y | zasadniczy: 1, 15,16; 31, 47, 78, 125, ... (%) 
TENGA 


sprzężony: 13, 27, 40, 67, 107, 174, 281, ... (X) 


By uniknąć, o ile to możliwe, otrzymywania wyrazów ujemnych w cią- 
gu złożonym, co mogłoby prowadzić do mimowolnych omyłek w obliczeniach, 
znalazłem następującą regułę empiryczną na tworzenie się jego wyrazów. 
Mianowicie: zgodnie z postacią ogólną (7) i zamieszczonym tamże ciągiem 
(7) — ciąg złożony (e) powstaje z ciągu (æ) przez spółezynniki ciągu (8), po- 
czynając od pierwszego; (£) powstaje z ciągu (8) przez spółczynniki ciągu 
(r), poczynając od trzeciego; ciąg (%) powstaje z (y) przez spółczynniki cią- 
gu (a), poczynając od pierwszego, i wreszcie (X) powstaje z ciągu (ò) przez 
spółezynniki ciągu ($), poczynając od trzeciego. Czyli: przy tworzeniu za- 
sadnieczych ciągów złożonych spółezynniki mnożącego ciągu zaczynają 
się od pierwszego; przy tworzeniu zaś sprzężonych złożonych, 
spółczynniki zaczynają się od trzeciego. 

Gdy będziemy rozważali postać algiebraiczną ciągów, przytoczony przy- 
kład przedstawi się jak poniżej 


; k, b kb kl, 2k4 31, «. (a) 
| 1—2k. 21—3k, 31—5k, 51 —8k, 81—13k, ... (B) 
> v, t, vt, v<-2t, 2v+- 3t, .. ©) (y) 

2U w, 24—3v, 3t—5v, 5t—8v, 8E—130, -... (2) 


*) Wiadomości Matematyczne. Tom XV. 1911. Str. 246—247, 
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(1—2wv)1—(21—3v)k, (24—3v31—(31—5v)k, (3t—5v)l- (54—8V)k,  ... — (8) 
| (v-+-t)(21—3k) —(v--26)(1—2k), (v+-24)(21—3k)—(2v-- 3t)(1—2R), 


S (9v4-31(21—3k)—(3v--50)(1—26) .. (0) 
= y i 
U kt—w, lU—(k--l)v, (k--I)t—(k--20pv, ... (x) 

| (31—5k)(21—3v) — (51—8k) (t--2v), (51—8k)(2t—3v)— 
—(81— 13k)(t—2v), (81 - 13k)(21—30)—(131—21k)(1—20), Q) 


Tym sposobem, w ciągach algiebraicznych (e), (£), (4) (X) mamy właś- 
nie tę konstrukcję wyrazów ciągów o określniku złożonym D=dyd,, która 
pozwoli nam wynaleźć elementy ciągów o określnikach d, i dy. 


Przypuśćmy teraz, że liczbę D kształtu 10m+1 mamy rozłożyć na 
czynniki pierwsze. Wymaga to przedewszystkim rozwiązania równania kwa- 
dratowego nieoznaczonego (5) — i wynalezienia początkowych wyrazów a i b 
ciągów, odpowiadających danemu określnikowi D 


b?—ab—a? =D (5) 


Podstawimy w ostatnim równaniu wyrażenie na przedwyraz U,=b—a 
i rozwiążemy przekształcone równanie (5a) 


U,b—a* =D lub inaczej U,(U, +a)=a? =D (5a) 


Wyrażenia na a i b przedstawimy w postaci funkcji zmiennych 0, i D; 
mianowicie PWDŁYEE: BARA 
m U +V 5 U,2—4D R 3U, +V50,—4D 

POLUS A ZHR an 2 


Znaki (—) przed pierwiastkiem stosują się do ciągu zasadniczego (licz- 
by mniejsze), znaki zaś -- do sprzężonego z nim. 

Chcąc z (10) otrzymać a i b wymierne, należy dobrać przez próbowa- 
nie takie całkowite U, abyśmy otrzymali pod pierwiastkiem kwadrat zu- 
pełny. , 

Dobór ten dokonywa się za pomocą tablicy kwadratów. Sądzę, że upro- 
ściłem go i udoskonaliłem do możliwych granie, a jednak czynność ta po- 
chłania 80%/, — 90%, czasu, potrzebnego do rozkładu liczby na czynniki 
pierwsze na podstawie mojej metody. 

Poszukiwanie wszystkich całkowitych pierwiastków na U, upraszeza się 
przez warunek 


(10) 


U,<VD CHI 
dowiedziony już uprzednio *). 


*) Wiadomości Matematyczne. Tom XV. Str. 241, 
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Znalazłszy to max.l/, na podstawie (11), obliczamy największą możliwą 
wartość wyrażenia podpierwiastkowego 


5(max. U„)=4D=Q (12) 


Podstawiając następnie liczby kolejne: max. Dę—1, max, Uy—2, i t. d., 
nie mnożymy już za każdym razem przez 5 kwadratu liczby kolejnej 
w tablicach, by odjąć potym 4D, lecz dobór próbny sprowadzamy do odej- 
mowania. Mianowicie, ponieważ różnica kwadratów dwuch liczb kolejnych 
równa się sumie tychże liezb, wystareza odjąć najprzód od znalezionego 
Q wielkość 5fmax. U,--(maa. Uę—1);, 


Q—5 fmax. U, (maa. U; — 1); = Q; (13) 


Każdy nowy odjemnik będzie mniejszy od poprzedniego o 10, gdyż 
suma dwuch liczb kolejnych różni się o 2 od sumy liczb sąsiednich. 

Odejmując tym sposobem liczby coraz mniejsze, prowadzimy tę czyn- 
ność tak długo, aż odejmnik okaże się większym od ostatniej reszty Q,. 
Wtedy należy sprawdzić działania arytmetyczne, posiłkując się bezpośred- 
nio zależnością 


5(min U,)?—4D= Qn (14) 


Gdy odejmowań jest bardzo wiele, należy to sprawdzenie stosować czę- 
ściej, np. w każdej zakończonej kolumnie pionowej. Wogóle ilość odejmowań 


U, 
będzie w przybliżeniu = n a + 50/9, czyli np. dla liczby D bliskiej miljo- 


nowi, należy wykonać około 105 odejmowań. 

Przekonawszy się ze sprawdzenia Q, według (14), że nie zaszły żadne 
pomyłki, porównywamy wszystkie reszty: Q, Qo Q» ... Qn z tablicą kwadra- 
tów, odznaczające te z nich, które się okażą kwadratami zupełnemi. Liczba: 
tych ostatnich powinna być równą 2. 

Ponieważ otrzymujemy stąd 2* różnych U, które czynią kwadratami 
zupełnemi wyrażenie podpierwiastkowe (10), zatym wszystkich ciągów typu 
Fibonacci'ego o tym samym określniku D będzie 2*+!, ze względu że każdemu VU, 
odpowiadają dwa ciągi sprzężone. Liczbę więc D można przedstawić w po- 
staci k--1 różnych form kwadratowych kształtu (5), czyli zawiera ona k--1 
różnych czynników pierwszych. 

Oezywistym jest, iż metoda niniejsza dostarcza zupełnego sprawdzenia 
działań arytmetycznych, wykazując jednocześnie liczbę różnych dzielników 
pierwszych w rozkładanej przez nas liczbie D. 

Gdy k=1, różnych U, będzie 2, i stąd dwa tylko dzielniki. Przy k=0, 
liezba D będzie pierwszą, gdyż 2%=1, czyli otrzymujemy jedyne U,, zamie- 
niające wyrażenie podpierwiastkowe (10) w kwadrat zupełny. z 

Jeśli się okaże, że żadna z reszt Q nie jest kwadratem zupełnym, ozna- 
cza to oczywiście, iż dana liczba D, jakkolwiek kształtu 1Om++1, nie może. 
być jednak określnikiem żadnego ciągu typu Fibonaccie'go. Zatym musi ona 
zawierać czynniki postaci 10m-+3. Taka liczba nieokreślnikowa (lub 
mieszana) będzie bezwątpienia złożoną, gdyż czynniki kształtu 10m--3 mu- 
szą znaleźć się eo najmniej w liezbie dwuch, by iloczyn ich otrzymał postać 
10m+1. 
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Znaleźliśmy zatym 2* różnych U,. Stosownie do ich wartości liczbowej, 
na zasadzie wzorów (10) układamy 2**! różnych ciągów typu Fib. Wystar- 
czy przytym obliczyć jedynie ich pierwsze wyrazy a, by rozwiązać zajmu- 
jące nas zagadnienie. Obieramy następnie dowolne dwie pary ciągów sprzę- 
żonych, bacząc by wyrazy a były liczbami możliwie małemi. Oznaczmy wiel- 
kości na a w jednej parze ciągów przez f i g, w drugiej przez h i ¿ i przy- 
równajmy do nich wyrażenia (9) na pierwsze wyrazy ciągów (e), (6), (%), 
(à). Otrzymamy wtedy cztery równania algiebraiczne 


(1—2v)1—(21—30v)k=f 
(v+V(2I—3k—(0-+-21)(1—2k) = g 


ki—lv=h (15) 
(31—5k)(2t—3v)—(51—8k)(t—2v)= i 
Rozwiązując je, znajdziemy 
p MADE vka 21 --39—3h—2i 
5 5 
| (16) 
nie eae 3f+79+3h+-2i 


Wszystkie liczniki w (16) dzielić się powinny przez 5. Łatwo spraw- 
dzić, że wszystkie one są porównalne względem modułu 5; zatym, gdy jeden 
okaże się wielokrotnością 5-iu, zachodzi to samo i dla reszty liczników. Wy- 
starczy więe najprostszy z nich (np. w wyrażeniu na kt) poddać podstawia- 
niu próbnemu, które w danym wypadku konieczne jest z tego względu, że 
oznaczenie literami /,g i h,i pierwszych wyrazów a w obranych przez nas 
ciągach — było w gruncie rzeczy tymczasowo dowolnym. Niewiadomo bo- 
wiem a priori z wyrażeń (9), które mianowicie z ciągów (a), (B), (y), (è) 
okażą się mnożącemi, a które mnożonemi, które zasadniczemi a które sprzę- 
żonemi, 

Jeśli licznik wyrażenia na kt nie dzieli się przez 5, to należy zastoso- 
wać litery f i g do drugiej pary ciągów obranych, litery zaś h i ¿ do pierw- 
szej. Gdy i to nie pomoże, należy w jednej lub obu parach ciągów przesta- 
wić te oznaczenia literowe z szeregu zasadniczego do sprzężonego i odwrot- 
nie. W rezultacie jedna z tych kombinacji okaże się niewątpliwie trafną. Co 
do tego łatwo się przekonać, wziąwszy zamiast licznika (/--4g--h-—t) po- 
równalne z nim względem . modułu 5 wyrażenie (f—g--h—t) i podstawiająe 
z (9) znaczenia na f, g, h, i. Otrzymamy wtedy ostatecznie 5v(k—l), czyli 
wielokrotność 5-iu. 

Dobór ten próbny wykonywa się bardzo szybko, gdyż nie obliczamy za 
każdym razem całkowitej wartości na kt, lecz orjentujemy się według ostat- 
nich cyfr liczb f, g, h, i, aby przekonać się, przy jakiej kombinacji oznaczeń 
ostatnia cyfra licznika kt będzie 5 lub 0. 

Znajdujemy tym sposobem na podstawie (16) wartości liczbowe na 
lv, kt, vk i tl. Będą one bardzo małemi liczbami w porównaniu z D, i bez 
trudności najmniejszej dobierzemy z nich natychmiast liczby całkowite na 
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k, l, v, t. będące czynnikami poprzednich. Że zaś są one jednocześnie wy- 
razami początkowemi ciągów (a) i (5) w (9), obliczywszy więc określniki 
d, id, tych ciągów, otrzymamy dwa gotowe czynniki liczby D. 

Gdy liczba D zawiera kilka dzielników pierwszych, wówczas jedna 
z liezb d) i d,, lub też obie mogą okazać się złożonemi. Z taką nową 
liczbą złożoną postępujemy zupełnie analogicznie, jak dotąd z liczbą D. Lecz 
otrzymane d, i d, są bez porównania mniejszemi od D, i na rozkład ich dal- 
szy wystarczy kilka minut czasu. 

Rezultat ten można osiągnąć inną drogą. Zamiast prowadzić rozkład 
liczb d, i d, za pomocą metody, zastosowanej do liczby D, zauważmy, iż ma- 
my do rozporządzenia k--1 ciągów, z których dotąd wyzyskaliśmy tylko 
cztery. Biorąc dwie inne pary ciągów, otrzymamy rozkład liczby D na dwa 
inne ezynniki d, i d,, przyczym naturalnie 


dyd,=dyd, . 


Znajdując największe wielokrotne czynników po obie strony tej rów- 
ności, będziemy już w posiadaniu 3-ch lub 4-ch dzielników różnych liczby 
B-4. Ło sd; 

Zdarzyć się może, że gdy wymnożymy wszystkie otrzymane czynniki 
proste rozkładu i podzielimy D przez ich iloczyn, to iloraz okaże się więk- 
szym od 1. Jest to zrozumiałe, gdyż potęga czynników, zawartych w D nie 
wpływa na liczbę odpowiadających mu ciągów. Zatym, iloraz ewentualny bę- . 
dzie się dzielił jedynie przez czynniki już znalezione. 

Cały przebieg czynności powyższej objaśnimy teraz na przykładzie licz- 
bowym. Rozłóżmy np. liczbę 0—638219 na czynniki pierwsze. 

Na podstawie nierówności (11), znajdujemy maa. UV,=798. Za pomocą 
wyrażenia (12) obliczamy resztę Q 


5 x 7982—4 x 638219=631144=0Q . 
Zgodnie z (13), należy odjąć od otrzymanego Q 
5(798--797)=7975 . 


Reszta będzie Q,=623169. Odejmujemy od niej 7965; od nowej reszty 
odejmujemy 7955 i t. d., jak to widać poniżej: 


631144 465769 304804 148249 
7975 7765 7555 7345 
623169 458004 297249 140904 
7965 71755 7545 7335 
615204 450249 289704 133569 
7955 7745 7535 7325 
607249 442504 282169 126244 
7945 7735 7525 7315 
599304 434769 274644 118929 
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599304 434769 274644, 118929 
7935 7725 7515 7305 
591369 427044 267129 111624 
7925 7715 7505 7295 
583444 419329 259624 104329 * 
7915 7705 7495 7285 
575529 411624 252129 97044 
7905 7695 7485 7275 
567624 403929 244644 89769 
7895 7685 7475 7265 
559729 396244 237169 * 82504 
7885 7675 7465 7255 
551844 388569 229704 75249 
7875 7665 7455 7245 
543969 380904 222249 68004 
7865 7655 7445 7235 
536104 373249 214804 60769 
7855 7645 7435 7225 
528249 365604 207369 53544 
7845 7635 7425 7215 
520404 357969 199944 46329 
7835 7625 7415 7205 
512569 350344 192529 39124 
7825 7815 7405 7195 
; 504744 342729 185124 31929 
7815 7605 190560 0 TIBS 
496929 335124 177729 24744 
7805 7595 7385 7175 
489124 327529 170344 17569 
7795 ŻA85.0 41 22385 7165 
481329 319944 162969 10404 « 
7785 7575 7365 7155 
473544 312369 155604 3249 * 
7775 7565 7355 
465769 304804 148249 


Wszystkich odejmowań wykonaliśmy tu 83. Zatym 
min.. U,=298—83=4215 . 


Otrzymane min, U, służy do sprawdzenia ostatniej reszty, t. j. powinna 
zachodzić równość 


5x 7152—4 x 638219—3249 , 
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` Przedtym jeszcze, w trakcie wykonywania odejmowań, sprawdzamy tą 
samą drogą reszty na końcu każdej kolumny. Że zaś mamy w kolumnie 21 
odejmowań, należy w wyrażenie podpierwiastkowe (10) wstawiać za U, 
liczby: 777, 756 i 735. 

Zwrócimy tu uwagę, że odejmowania powyższe możemy zastąpić przez 
dodawania, wychodząc z min. (05=715, które łatwo wyszukamy przy pomocy 
tablicy kwadratów, i orjentując się według wskazanej powyżej ogólnej przy- 
bliżonej liczby wszystkich odejmowań oezekiwanych. Pierwsza dodajna będzie 
się równała 5(715--716)=7155; każda następna będzie większa o 10 od po- 
przedniej. Szereg kolejnych dodawań powinien nas doprowadzić do sumy 
Q=631144. 

Po sprawdzeniu odejmowań (czy też dodawań), przekonywamy się, że 
śród naszych reszt istnieją cztery kwadraty zupełne. Oznaczone są one gwiazd- 
kami i równają się: 487%, 323*, 102? i 57%. Odpowiadające im przedwyrazy 
U, będą: 747, 729, 716 i 715 — eo łatwo się wynajduje przez obliczanie 
ogólnej liczby odejmowań, wykonanych od samego początku. 

W danym przypadku 4=2?, t. j. k=2, a więc liczba D zawiera trzy 
czynniki pierwsze. Wszystkich ciągów o określniku D=638219 będzie 8; otrzy- 
mamy je łatwo na podstawie wzorów (10). Mianowicie 


i f 130,5 877, ... f 203, 932, ... 
HTI 1 614, 1364, .. GRY | 526, 1255, ... 17 
ię | 307: 1025, -. ąz | 32% 1084, 

6 | 409, 1125, ... 8 aęk, 380 SIŚD a 


Z pomiędzy (17) obierzmy dwie pary górne ciągów, i zgodnie z układem 
równań (15)— oznaczmy: f=130, g=617, h=203, 1=326. Przy pomocy wzo- 
rów (16) znajdujemy 
; v=252; kt=455; vk=90;, tl=1274. 

Stąd oczywiście 


k==5, v=18, t=9l, 1=14 
czyli poszukiwane ciągi czynnikowe (œ) i (y) w (9) będą 
SUPER d4= 101 
18, 91, 109, 200, ... d+==6319. 


A więc 101 jest jednym z czynników pierwszych liczby 638219; iloczyn 
dwuch innych równa się 6319. 

Obierzmy z pomiędzy (17) dwie pary ciągów lewych, i oznaczmy: f=307, 
g=409, h=617, i=130. Wzory (16) dadzą 


lv=—131; kt=486; vk=—54: tl=1179. 
Skąd 
v=—]; k=54; t=9; l=131 
a same ciągi czynnikowe będą 
54,43); 185, "316, ... d;=7171 


1. 00 8, Ww 04 BU 
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Otrzymaliśmy więc drugi czynnik poszukiwany 89. Trzeci znajdzie się 
już łatwo jako 71. Dość jest podzielić d, przez d, lub też dą przez d,. 

Gdybyśmy zaś rozkładali liczbę d =6319 za pomocą metody psi to. 
rezultat byłby równie szybki. 

Odpowiednie jej max. (/,=79 i reszta `Q 


| 5x 79%—4x6319=5929=Q. 
Wykonamy następnie odejmowania, poczynając od 5(79--78)=785 


5929 * 2849 
785 745 
5144 2104 
775. 735 
4369 1369 * 
765 725 
3604 644 
755 
2849 


Dwie. z reszt powyższych są kwadratami zupełnemi: 5929 = 77? 
i 1369=37?. Odpowiednie U, równają się: 79 i 73 i na podstawie (10) ciągi 
odnośne 
f 1, 80, 81, .. 18, 91, 109, 


l 78, 157, 235, 55, 128, 183, 
Oznaczając: f=1,.9=78, h=55, 4=18, znajdujemy według (16) 
o =A REE PKE708 2150 


U,=79 0,=23 | 


czyli 
k=7, v=l,. t=10, I=15 


i napiszemy ciągi czynnikowe, obliczając ich określniki 
T ASSESI y A d 71 
1,210, 10ZŁ E dą =89. 


Wybór tego lub innego sposobu zależy od wielkości liczby rozkładanej, 
od ilości ogólnej zawartych w niej czynników pierwszych, większej lub też 
mniejszej złożoności pierwszej pary d, i d, znalezionych dzielników. Wogóle 
łatwo zorjentować się, który z nich należy zastosować. 

Z przytoczonego przykładu liczbowego widoczne są korzyści zastosowa- 
nia metody tu wyłożonej. Wszystkie czynniki pierwsze otrzymujemy jako 
liczby gotowe, nie zaś za pomocą dzielenia próbnego liczby danej przez 
liczby, znalezione drogą uboezną, jako możliwe jej czynniki. Ilość zawar- 
tych w liczbie badanej czynników pierwszych wskazujemy a priori. Wreszcie 
co do szybkości rozwiązania, to przeciętnie robi się cztery odejmowania 
na minutę. Posiłkując się zaś arytmometrem, można to przyspieszyć 2 do 
3-ch razy, przyczym jedną osoba działa arytmometrem, a druga porównywa 
otrzymane reszty z tablicą kwadratów. Gdy liczba okaże się złożoną, to na 
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wszelkie czynności pozostałe, jakoto: ułożenie odpowiadających jej ciągów, 
układu równań (16), znalezienie ciągów czynnikowych, dość jest 10 — 15 mi- 
nut, zależnie od wielkości liczby rozkładanej i liczby jej czynników pierw- 
szych. Np. w celu zbadania liczby Czebyszowa D=8520191 należy wykonać 
307 odejmowań, co wymaga. około 80 minut, (z arytmometrem najwyżej 
30 min.). Okaże się przytym, że jedna tylko reszta 58081 jest kwadratem 
zupełnym, więc liczba badana jest pierwszą, i żadnych więcej już z nią 
ezynności niema potrzeby wykonywać. Również np. dla liczby 1 111111 znaj- 
duję w ciągu 30 min. bez arytmometru czynniki jej 239 x 4649. 

Największą zbadaną przezemnie liczbą był wyraz 37-y szeregu 1, 3, 4, 
7, 11,... który proponuję nazywać szeregiem R, by go wyodrębnić ze wzglę- 
du na pewne specjalne własności. 


Uk — UF + UE = 54018521 . 


Liczba ta jest pierwszą. Wszystkich odejmowań było 772, eo wykona- 
łem w ciągu 70 minut przy pomocy arytmometru. Tablica kwadratów 
I. Claudela pozwala na rozkładanie liczb określnikowych do wielkości 
około 100 miljonów, i eała metoda moja zależną jest od dokładnej a obszer- 
nej tablicy kwadratów. | 

Pozwolę tu powołać się na taką powagę w teorji liczb, jak E. Lucas 
(Théorie des Nombres. 1891. Str. 36), który powiada: „Tablica kwadratów 
pozwala na ogromne uproszczenie rachunków arytmetyki; nie wątpię, że ta- 
blica ta była w rękach Fermata cudownym narzędziem obserwacyjnym 
w jego poszukiwaniach arytmetycznych i punktem wyjścia jakiejś szybkiej, 
dziś prawie nieznanej metody rozkładania liczb na czynniki pierwsze“. 

Otóż zdaje mi się, iż rezultaty metody mojej wskrzeszają w tym wzglę- 
dzie powagę tablicy kwadratów i uchylają rąbek owej tajemnicze) zasłony. 

Do dalszej pracy nad tą kwestją wzywam matematyków kompetentniej- 
szych odemnie, sam badań tych nie zarzucając. Jest do zrobienia jeszcze 
bardzo wiele; mianowicie należy: 1) wynaleźć formę ogólną dla liczb nieo- 
kreślnikowych, poczym łatwo będzie stworzyć metodę ich rozkładu; 2) zna- 
leźć krótszy od mego sposób oddzielania czynników określnikowych od nieo- 
kreślnikowych w liezbie badanej; 3) może uda się komu odkryć jeszcze szyb- 
szy sposób na wynajdywanie całkowitych U, w wiadomym nam wyrażeniu 
podpierwiastkowym, by uniknąć ogromnej ilości odejmowań, koniecznych, jak 
dotąd, przy rozkładzie liczb bardzo wielkich. 

R. Niewiadomski. 


www.rcin.org.pl 


